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The primary problem in the calculation of the differential 
scattering cross section for hydrogen molecular ions is the derivation 
of an approximate expression for the probability that the hydrogen 
molecular ion will be oriented, prior to impact, in such a way as to 
avoid disruption during scattering. Using a simple model, an approxi- 
mate expression for the orientation probability is obtained. 
Differential cross sections for the scattering of hydrogen molecular 
ions from carbon and krypton are calculated using classical scattering 
theory and the derived orientation probability. Impact energies of 
100 to 1000 electron volts are considered. In addition differential 
cross sections for protons scattered from carbon and krypton are 
presented. 
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1 . INTRODUCTION 



Bates, Cook, and Smith, have presented a classical theory of 
ion-molecular rearrangement collisions at high impact energies^^ . 
Using proton-methane collisions as an example a brief outline of the 
theory follows. 

In Figure 1 m^ is the impacting proton, mg is one of the 
hydrogen atoms of the methane molecule and is the carbon atom of 

the methane molecule. In this collision nu = mg « M^. For capture 

/ // 

to occur re V g • From these two conditions and the conservation of 
of linear momentum and energy it can be shown that 0g re 90° and 
0^ re 45°. Thus, for this case, capture of the hydrogen atom of the 
methane by the incident proton to form hydrogen molecular ions may 
occur at an angle of approximately 45°. Bates, Cook and Smith have 
calculated a cross section for the production of hydrogen molecular 
ions by this process^"^. 

An experimental difficulty may arise in a measurement of the 
cross section for the production of Hg by this double scattering 
process if the incident proton beam contains some significant number 
of Hg . If the cross section for scattering of Hg from the carbon 
atom is much larger than the cross section for Hg production by capture, 
then a measurement of the production cross section would be impos- 
sible unless a pure source of protons were used. In order to 
determine the relative magnitudes of these two differential scatter- 
ing cross sections a calculation of the differential cross section 
of Hg ions scattered from heavy atoms is presented. 
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Figure 1. Schematic illustration of the capture process 



2. THE MODEL 



The hydrogen molecular ion is linear with an inter atomic spacing 
of approximately two angstroms. For randomly oriented impacting 
ions this implies a possible difference in impact parameter for the 
individual atoms of as much as two angstroms. If the difference in 
impact parameters is great enough, the velocities of the two atoms 
after scattering may differ by enough to cause their relative energy 
to exceed the binding energy of the molecule. The molecule will 
then be disrupted and can not be considered an H 2 scattering event. 

Even at zero difference in impact parameter disruption of the 
molecule may occur because of the following consideration. In a 
spherically symmetric potential, such as that used in this calcu- 
lation, each particle remains in a plane, defined by the scattering 
axis and a line from the scattering axis to the particle, during 
and following the collision. If the angle between the planes is 
great enough the difference in the velocities of the two atoms, 
after scattering, may be such that the molecule is disrupted. 

Figure 2 shows the relationship between the scattering planes 
and the scattering axis. £4> is the angle between the scattering 
planes. 

The condition that the molecule remains intact may be stated as 
D > % /i(£V ) 2 where D is the binding energy, £V is the relative 
velocity of the two particles and % /i(£V ) 2 gives the relative energy 
of the two particles. Using the symbol AV for the maximum allowable 
£V the limiting condition for nondisruption becomes: 

D = k H (AV ) 2 . 
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Figure 2. Illustration of scattering model. 



Thus both the difference in impact parameter Qp and the angle between 
the scattering planes must be such that the above condition is 
satisfied. This limitation on the orientation of the impacting 
molecule leads to an orientation probability that must be considered 
in the calculation of the differential scattering cross section for 
hydrogen molecular ions. 

For high impact energies and small impact parameters each atom 
of the may be treated as a proton. That is, the effective charge 
may be taken as one for each particle. The calculation of the differ- 
ential cross section for Hj£, ions may then be approached in the follow- 
ing way. For a given energy and impact parameter calculate the 
scattering angle, and the differential cross section for that angle, 
for one proton whose velocity is V 0 , the velocity of the h£ . Then 
determine the probability that the "second" particle is oriented with 
respect to the "first" particle in such a way that its scattering 
angle is within A 0 of 8 and its scattering plane makes an angle of 
less than A $ with the scattering plane of the first particle. The 
differential cross-section for Hg scattering will then be given by: 
O h *(B) = CT H t(0) where P is the orientation probability. This is 
then done over a range of impact parameters and for several energies. 

3. CALCULATION of a Ht (0) 

The calculation of a H< (0) for the "first" particle was done using 
the method suggested by Bates, Cook and SmitlJ"^. The potential used 
for the calculation was a shielded potential of the form, 

V(R) = ^ (Me*^ Ne' + Le'*) due to Byatt^^. Calculations were made 
for both carbon and krypton. Byatt lists several potentials in 
reference 2, all of which may be conveniently used in the computer 
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program in appendix 2. The cross section was calculated for impact 
parameters from 0 - 5 a 0 in increments of 0.06a o (a 0 is the radius of 
the Bohr orbit) and for energies from lOOev to lOOOev in steps of 
lOOev. The details of the calculation are shown inappendix 1. 

4, CALCULATION OF THE ORIENTATION PROBABILITY. 

The condition AV =/ 4D may, for convenience, be put in terms 

\%/ E^ 

of the energy of the . E^ = or ^ = -S 2 . then 

V Q 

AV = ^4D j ^ V 0 . Th’is limiting condition implies a relationship 

between QQ and £<i> such that AV= F(£0, £0 is related to £p by 

£0 = p£p/a„<(0) SIN0 in the center of mass coordinate system. 

Therefore AV = G(£<l>,Gp) gives a relationship between £<£ and £p. Since 
AV is a constant for a given energy and impact parameter, there exists 
some relationship between Qp and such that £4> = f(£p). £<£> = f (£p) 

defines a curve in a plane perpendicular to the scattering axis. The 

. 1 j 

second particle must pass through the area on the plane bounded by 
this curve if disruption is to be avoided. 

The randomly oriented can be considered as defining a sphere 
the center of which is the "first 11 atom and whose radius is JL 9 the 
inter atomic spacing. 

In Figure 3, A is the area on the surface of the hemisphere con- 
tained within the surface defined by Qp = f (£40 • The ratio of this 

j 

area to the area of a hemisphere of radius & gives the fraction of 
the ions which are oriented in such a way as to avoid disruptions. 
Thus the orientation probability, P is given by: P = A/27 TJL 2 . No exact 

solution has been obtained for the orientation probability. 

As an approximation to the orientation probability the following 
simplifications were made. 
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Figure 3, Illustration of orientation probability model 



Assuming the target atom to be infinitely massive, there will be 



no change in the magnitude of the velocity of either atom of the 
hydrogen molecular ion. Since they then have the same initial speed 
V 0 , Vl = V 2 = V 0 . When £<i> = 0 the particles are scattered in the same 
plane as shown in Figure 4a. Then, assuming 40 is small 
AV « Vo 40 or 40 . When £0 = 0 (wfricfr implies £p = 0) 

the relationship between the individual atoms is as shown in Figure 4b. 



Then AViS (V. SUB) or ^ ‘ 

"tr 



* 






these approximations 



then give independent maxima for £<£ and £0. £0 is related to £p by 

£p = y(9) SIN0 ^0 or Zip = cr H+ (0) SIN0/ 4D 



ti) 



The above restriction on 4p implies that the second patticle must 
be contained between two coaxial cylinders whose axis is the scattering 
axis and whose radii are p - 4p and p + 4p respectively. The condition 




puts a limit on the angular separation of the scattering planes. The 
limits on 4*£ and 4p then confine the second particle to lie on the 
surface of a hemisphere of radius l contained between the two scattering 
planes and the two cylindrical surfaces. Figure 5 illustrates. 

This area, on the surface of the hemisphere, divided by the area 
of the hemisphere then gives an approximation to the orientation 
probability. 

As a further approximation, this area was assumed small and taken 

to be the plane area given by A « p(2A<l>) (2 4p) ; which then gives, as 

an estimate of the orientation probability* P « 2pAjjE£p 4 

Til 
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Figure 4. Illustration of the approximation to the orientation probability 



I 

t' 

I 





2T 



18 



Figure 5. Illustration of the approximate probability area 



5. RESULTS AND DISCUSSION 



Everhart, et al, have calculated the differential cross section 

[3] 

for scattering from a screened Coulomb potential . The results are 
presented in terms of a, the screening length, given by: 

a = [Zi* + z? J % 

and b, the collision diameter, given by: 

b = a Zp e 3 

hVr 

where \i = the reduced mass of the system. 

V = the relative velocity of projectile and target. 

As a check on the calculation of the cross section for proton 
scattering the computer program used for the calculation was run using 
the screened Coulomb potential of carbon and a proton impact energy 
such that the ration b/a was equal to 1. The results of this calcula- 
tion were within 2.34% of those of Everhart et al. 

The actual calculations were made using the Byatt potentials of 
carbon and krypton with proton impact energys from 50 to 500 electron 
volts. The results are presented in the graphs of Figures 6 and 7. 

The cross section for the scattering of the hydrogen molecular 
ion for impact energies of 100 to 1000 electrons volts (corresponding 
to proton energies of 50 to 500 electron volts) then is the cross 
section for proton scattering multiplied by the orientation probability. 
These cross sections are presented in the graphs of Figures 8 and 0. 

Under the assumption that the speeds of the two particles are the 
same following scattering, the LQ used in the orientation probability 
approximation is due only to the difference in the directions of the 
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Scattering angle-degrees (lab. coords.) 

Figure 6. Graph of the logarithm of the differential scattering cross 
section versus scattering angle for protons on carbon. 
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Figure 7. Graph of the logarithm of the differential scattering cross 
section versus scattering angle for protons on krypton. 
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Figure 8. Graph of the logarithm of the differential scattering cross 
section versus scattering angle for hydrogen molecular ions on carbon. 
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Figure 9. Graph of the logarithm of the differential scattering cross 
section versus scattering angle for hydrogen molecular ions on krypton. 
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velocity vectors. The speed of a particle scattered through an angle 
0 is given by: 



V = V. 



i- 2M A 

(V M if 



(1- COS0) 



where 



V = the speed of a projectile particle scattered through 
an angle 0 in the C.M. Coordinate system, 

M^= Mass of the projectile particle, 

M fc = Mass of the target particle 



and 



V 0 *f initial speed of the projectile particle. 

If the speeds of particles scattered through angles 0 and 0 + C0 
differ by enough that GVg, resulting from this difference, is greater 
than A V, the maximum for non disruption, then the Hg will be disrupted 
due to the speed difference alone. If 40 s »40, then £0 (due to the 
difference in direction only) will be the more stringent condition and 
the error introduced by ignoring the difference in speeds may be assumed 
small. 
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2M M, 

The factor t n , for scattering from carbon, is approximately 
(M H +M t ) y 

. The factor (1- cos0) is always less than one for scattering 



angles between zero and 90 degrees. Thus the expression for V above 
can be expanded in a series. Taking only the first two terms gives. 



V = V„ 



M t Mh 



then £V = -V 0 



M 



t “h 



K + \ r 

sin© C0 



(1- COS0) 



’] 



when £V 5 = AV and £0 5 = A& then |^0 C | 

VtsiveV'l 
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but A V = 



_/4D 




1“h 



and V 0 = 



■m 




, which is just the £Q due to the differ - 



Therefore A V = /4D 

V ° (%■ 

ence in direction t>f the velocity vectors. 

(K+m \ 2 

Therefore £0 = V— - — . For carbon, then Z0^ ^ 14 Z^0. 

M^M^sinQ 

This implies that the will be disrupted due to the difference in the 
directions of the velocities of the two particles at a £0 which is much 
less than the £0^ required to cause a disruption due to a speed differ- 
ence only. 

Because there is some difference in the speeds of the two particles 
the relative energy of the particles is somewhat greater than that 
predicted with the speeds assumed equal. Thus disruption will occur 
more often than predicted by the model used for the calculation. 
Therefore, the assumption of equal speeds tends to make the probabil- 
ity factor somewhat too large. 

The use of separate maxima for Z0 and A<£ introduces an error in 
that the final probability expression contains the product of the 
maximum values, each of which was derived under the assumption that 
the other was zero. This error tends to make the probability factpr 
too large. 

In view of the small values of P, (representative values are 
listed in Table 1.) the assumption that the area on the surface of the 
hemisphere can be approximated by a plane area would seem to be 
reasonable for, at least, the higher angles and energies. 

The conclusion is that the probability factor is in general some- 
what too large thus causing the calculated cross sections for 
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Energy 
(in ev) 


0 (Lab) 


Carbon 

P 


Krypton 

6 (Lab .) P 


100 


83.65 


1.95xl0 _J 


87.88 


5.08x10 


100 


45.62 


7.72xl0“ 3 


45.23 


2.22x10 


100 


10.68 


2 . lOxlO -1 


10.32 


6.38x10 


500 


88.88 


9 . 18xl0~ 5 


83.45 


2.68x10 


:'0 500 


42.35 


7.08xl0 -4 


45.39 


1.66x10 


500 


10.37 


2. 30xl0 -2 


10.10 


7.00x10 


1000 


75.20 


1.72xl0 -5 


86.47 


6.04x10 


1000 


47.62 


7 . 12xl0~ 5 


46.80 


3.12x10 


1000 


10.35 


7.46xl0 -3 


10.13 


2.32x10 



Table 1. Representative values of the orientation probability. 
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2 

1 
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2 
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hydrogen molecular ion scattering to be too large. 

In Figure 10 the curve marked CT h *( 45° ) is the differential scatter- 
ing cross section for from carbon at a scattering angle of 45° as 
a result of this calculation. The impact energies plotted are from 

lOOev to lOOOev. The curve labeled aC^+Cl^ > CHq+H^) is the cross 

section for the capture of a hydrogen atom from the methane molecule 
by an impacting proton. This curve is a result of calculations based 
on the reaction cross section as derived by Bates, Cook and Smith^^. 
Figure 10 shows the cross section for scattering of from carbon at 
45 g to be somewhat more than an order of magnitude greater than the 
cross section for the reaction H + 4-CI^ — >Ctfe+H^. As a result, in 
order to measure the reaction cross section it would appear that a 
proton beam free of will be necessary. 
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LOQq Energy (in electron volts) 

Figure 10. Comparison of reaction and scattering cross 
sections . 
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APPENDIX I 



CALCULATION OF (©) 



The classical differential scattering cross section in the center 
of mass coordinate system is given by: 

crce) = <? (l) 

where Pis the impact parameter and B is the scattering angle in the 
center of mass coordinate system, ©is related to(?by: 0= 7T- 2 €U (2) 



where 



u = 



r 



CO 



A 



AT 






2 *) '/2 



in which is the largest positive root of 

I _ V(tL _ - n 

I yzMV} r* " u 

and Ai - the reduced mass of the system 
V(f) = the scattering potential. 



(3) 



(4) 



V$ = the velocity of the impacting particle in the lab 

coordinate system when the target particles are considered 
to be at rest. 

[2] 

The potential as given by Byatt, is 

V(r) = (oe ar+ Ee tr + Fe cr ) 

after converting V(f) to electron volts and letting fi(( 0 there results: 

V(R) =■ ^■(Ae° ,R + se pR + ce™) 





For Carbon 


For Krypton 


A= 


204.136 


328.635 


B= 


-71.856 


588.600 


C= 


31.029 


63 . 765 


a= 


-1.700 


-1.086 


fi- 


-2.895 


-4.961 


y= 


-8.808 


-26.110 


The 


quantity can be put in terms of the 


energy of the H* 



(5) 

( 6 ) 



in electron volts in the lab. coordinate system. 




2 

M«| 



2 Eh? . Eft 
2 
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in which the mass, N1 r ,is roass of the target atom. 
Mh + Alf * 

Define ET = L ( - & .. $ *-) £mL 
' 2 - \ Wm+Mt / 



( 7 ) 



Define V = 


ma. 

ET 




(8) 


Then fi 0 is the root of 


v. - n 

r /?* “ ° 






1 " 


(9) 


The root of 


this equation was found using the Newton- Raphson method 




which is 


fl.t, = ft-* 


_ fc«.) 


(10) 


where 


m = i 


y. _ <** 

/? fi* 


(ID 


and 


f(R) = v 


APV. t ££* 


(12) 


in which DV 


is defined by 

DV = 

dR 


&Pe f *+ cre ,f ) 


(13) 


then 


ET 




= »„ 


_ ffi* -Vfr. - 

(V - /?*DV +2<py/?,) 


(14) 



The integral, U, was evaluated using the Gauss-Mahler Quadrature 

[i] 



formula as employed by Bates, Cook and Smith . The result is 

^o5^)^(co^) +(cosf)t(cos$) +Ccojg)^(co5ff)J 

where - 4 a 1 “ 






in which the 



Xi s CO sfg ) Cos $ Avd 

To calculate jsL^-j differentiate equation 2 with respect to 



(15) 

(16) 

(17) 



Id© 

Ide 




(18) 
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where 



dj£ - 2k + 
de 9 Rode 



(19) 



_ 

3$ 



“ A ♦(< w f ) f£ (c olf) *(co J ff coi 



in which ^ _ , v - , - 1“ 

§f * «.ex?[fi* -va -rttf] 

and the X^, are as defined previously. 

dR 0 

is obtained by differentiating 
the result is 

dRo _ 2 . e 



de 2e 2 //? o -R C DV + V 

& - a [^®ir s5) ♦<*««&'•*#>] - % 



in which 



= ^ ' v/f ” x ‘ ' *’** ] [x " % - ^ ? ] 



where the are as defined previously, 



( 20 ) 



( 21 ) 



( 22 ) 

(23) 

(24) 
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FORTRAN 60 PROGRAM FOR THE CALCULATION OF THE DIFFERENTIAL 
SCATTERING CROSS SECTION FOR HYDROGEN MOLECULAR IONS FROM 
HEAVY ATOMS 
PROGRAM VELVA 
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EXPLANATION OF PROGRAM 



STATEMENT INDEX 



1-5 


Places and prints column headings for print out. 


6-8 


Dimension statements. Self explanatory. 


9-23 


Do loop for graph title. 


24-29 


Constants for the potential used. See equation 
No. 6 , Appendix 1. May take on any values in 
accordance with potential used. 


30 


Sets the initial value for the number of 
iterations. See statement index 74 and 75. 


31 


Sets the initial value of the impact parameter 
such that RH0(1) = Q. See statement index 60 


32 


Sets the initial value of theta, the scattering 
angle in C. M. Coordinates. 


33 


Sets the initial value of the energy of the 
hydrogen molecular ion in the LAB. Coordinate 
system, EV, such that EV(1) = lOOOev. See 
statement index 37. 


34, 39, 40, 41, 
42, 43, 44, 45 


Give the draw program the proper parameter to 
draw 10 curves on one graph. 


35 


Begins outer do loop for 10 different energies 
defined by statement index 33 and 37. 


36 


Subscript relationship to achieve desired energies 


37 


Energy defining relationship. May be varied for 
desired energies. 


38 


Sets the initial value of R for the iteration 
to find the root of equation 4, Appendix 2. May 
be based on b, the collision diameter. (Pag el 9 ) 


47-57 


Label program for the labeling of curves with 
appropriate energy. 


58 


Begins inner do loop for 101 values of the impact 
parameter, RH0, as defined by statement index 31 
and 60. RH0 is in units of a° 


59 


Subscript relationship to achieve desired values 
of RH0. 
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60 

61-63 

64 

65 

66 

67 

68 

69 

70 

71 

72-73 

74 



Defining relationship for RHO. May be varied 
as desired. 

Evaluates the exponential terms of the potential 
used. 

Mass of the projectile particle in a, m. u. 

Mass of the target particle in a, m, u. 

Evaluates the relative energy factor of equation 
4, Appendix 1 in terms of EV. 

Evaluates V as defined by equation 8, Appendix 1, 

Evaluates equation 11, Appendix 1, for the 
iteration. 

Evaluates DV as defined by equation 13, Appendix 
1, for the iteration. 

Evaluates equation 12, Appendix 1, 

Computes a new value of R for the iteration* 

Compares the last two values of R, 

Counts the number of iterations. 



75 



Stops the program if finding R requires over 
10000 iterations. 



76 Replaces R with the last value computed by 
statement index 71. 

77 Repeats the iteration if statement index 73 has 
not been satisfied. 

78 Continues the program if statement index 73 has 
been satisfied. 

79-81 Insures that the proper value of R has been used 

to compute V and DV. 

82 Stops the program if R is less than RHO and exits 

the inner do loop if R = RHO. 

83-85 Evaluate the cosine functions of equation 15, 

Appendix 1. 

86-100 Evaluate the g (X.) as defined by equation 16, 

Appendix 1. 
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101 

102 

103-105 

106 

107 

108-113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124-127 

128-132 



Evaluates the integral U as defined by equation 
15, Appendix 1. 

Computes theta the scattering angle in the C. M. 
coordinate system for the given energy and impact 
parameter by equation 2, Appendix 1. 

Computes the 3 p as defined by equation 21, 
Appendix 1. 

Evaluates 30/9* as defined by equation 20, 
Appendix 

Computes 
Appendix 

Computes 
Appendix 

Computes 
Appendix 

Computes 

1 . 

Computes 

1. 

Computes 

scattering of protons at the previously computed 
angle by equation 1, Appendix 1. 

Computes A pas given on page 16. 

Prevents dividing by zero in statement index 122 
when theta is 180. 

A dummy statement to give a value of A $ for 
theta = 180. 

o 

Skips statement index 122 when theta = 180. 

Evaluates as given on page 16. 

Evaluates^! as given on page 16. 

Compute the percent difference in velocities 
after scattering as described on page 24. 

Convert the scattering angle, theta, to LAB 
coordinates and degrees. 



1 . 

d & as defined by equation 22, 

1 . 

the as defined by equations 24, 

1 . 

SU /Sfto as defined by equation 23 

1 . 

dU/dp as defined by equation 19, Appendix 
!<■•/*. I as defined by equation 18 r Appendix 
the differential cross section for 
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133 


Gives the inter atomic spacing 1 . 


134 


Computes the value of the orientation probability 
as defined on page 


135-136 


Limits the value of P to 1. 


137 


Computes the differential cross section for 
H 2 scattering as given on page 13 # 


138 


Gives the proton energy for printing when that 
is of interest. 


139 


Prevents the log. function from having a negative 
or zero argument in statement index 140. 


140 


Computes the logrithum of the cross section for 
the draw program. 


141-142 


Limits the values of the scattering angle printed 
to those between 2° and 90° . 


143 


Prints results for each value of RHO assigned. 


144-150 


Self explanatory. 
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